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Prediction for the transverse momentum distribution 
of Drell-Yan dileptons at PANDA 

0. LinnykQ K. Gallmeister, S. Leupold, and U. Mosel 

Institut fur Theoretische Physik, Universitat Giessen, Germany 
(Dated: February 2, 2008) 

We predict the triple differential cross section of the Drell-Yan process pp — > X in the 
kinematical regimes relevant for the upcoming PANDA experiment, using a model that accounts 
for quark virtuality as well as primordial transverse momentum. We find a cross section magnitude 
of up to 10 nb in the low mass region. A measurement with 10% accuracy is desirable in order 
to constrain the partonic transverse momentum dispersion and the spectral function width within 
±50 MeV and to study their evolution with M and \fs. 
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Using the model described and tested in an earlier pub- 
lication [J, we calculate the p^-spectrum of the Drell- 
Yan pairs in pp collision at the centre of mass energy 
= 5.5 GeV, which is an interesting, not yet experi- 
mentally explored kinematical regime. The planned ex- 
periment PANDA 0] at future GSI facility is to probe 
the Drell-Yan process at this energy. 

The conventional perturbative QCD (pQCD) approach 
to the calculation of hard scattering cross sections as- 
sumes that the partons are collinear and on-shell, i.e. 
both the virtuality and transverse momenta of the quarks 
and gluons inside hadrons are neglected and the partonic 
four-momentum is replaced with it's light cone projection 
(fc + ,0, 0). This approximation to the kinematics does 
not effect the cross section of fully inclusive processes 
like DIS, but it leads to unrealistic results for the cross 
sections of more exclusive processes. Intrinsic transverse 
momentum of quarks has been proven to be important 
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FIG. 1: Relevant initial state interactions that can build up 
intrinsic transverse momentum and width of quarks in the 
proton in the Drell-Yan process. 
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for the interpretation of various cross sections and asym- 
metries 3]. The virtuality of quarks and gluons in the 
proton has attracted considerable attention recently. As 
shown in pj, by taking into account the virtuality dis- 
tribution of partons one can explain the K-factor type 
corrections to the transverse momentum distribution of 
Drell-Yan pairs. The double differential cross section of 
the Drell-Yan process da/dMdxp can be well described 
by pQCD in the conventional collinear on-shell approach. 
However, the pt distribution, i.e. the triple differential 
cross section da / dMdxpdpT in every order of collinear 
pQCD is unphysically sharp. One has either to resum all 
orders by the Sudakov method, introducing additional 
rcgularization parameters, or to correct the kinematics 
by taking into account quark transverse momentum and 
virtuality. The latter procedure leads to a very good de- 
scription of data already in the leading order of pQCD 
for the elementary process [IJ. The importance of ac- 
counting for the virtuality of partons in other high en- 
ergy processes was stressed, for instance, by O. Benhar 
and V. R. Pandharipande |4| and, recently, by J. Collins 
and H. Jung 

The generalization of the usual parton distribution 
function, which depends on the virtuality and kx of the 
quark in addition to it's longitudinal momentum k + , is 
sometimes called "double unintegrated parton density" 
or "parton correlation function". But the terminology is 
still not fixed, and one of the reasons is that the gauge in- 
variant definition of these distributions in terms of parton 
field operators is not known yet. It is an open question, 
whether factorization still holds in the generalized case, 
in which the soft part of the cross section depends on 
the full quark momentum {k + , k~, hp)- Also, the evolu- 
tion of these parton distributions has not been studied 
so far. Before the factorization theorem is proven, one 
can calculate measurable cross sections using a model 
that assumes generalized factorization. Of course, the 
predictive power of this model has to be tested. Com- 
parison of such calculations to the data in allowed 
us to fix parameters of a simple factorized parametriza- 
tion for the unintegrated quark density in the proton. In 
the present paper, we use the same parameters to cal- 
culate the cross section in the kinematics of the future 
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PANDA experiment in order to obtain a first estimate 
for the counting rates and p^-distributions expected at 
y/s = 5.5 GeV, M = 1 — 4 GeV, which is on the edge of 
the non-perturbative regime. To this aim we also com- 
pare with the PYTHIA event generator [6j. 

In [l| , a model to calculate the initial state interaction 
(ISI) effects in the Drell-Yan process was proposed. This 
model assumes, firstly, the generalized factorization and, 
secondly, that the soft part of the cross section can be 
approximated by a product of functions of the quark k + 
momentum, the transverse momentum, and the virtual- 
ly m 2 = k + k~ — k\. We use a Gaussian fc^-distribution 
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so that the mean squared partonic intrinsic transverse 
momentum (k\) = AD 2 . 

The authors of |4j pointed out the analogy between 
the more conventional collinear approach and the plane 
wave impulse approximation (PWIA) of many body the- 
ory. Corrections to the PWIA are due to initial and 
final state (FSI) interactions. ISI of active and specta- 
tor quarks generate quark virtuality and intrinsic trans- 
verse momentum (see Fig. ^ . Calculations in the quark- 
diquark model show that the partonic ISI corrections 
to the cross section of the unpolarized Drell-Yan pro- 
cess are not small 7]. However, one has to go beyond 
the single gluon exchange approximation of in order 
to study the ISI effects quantitatively. At present, the 
quark virtuality distribution cannot be calculated from 
first principles QCD. Using analogy to many body the- 
ory, we parametrize the quark virtuality distribution as 
a Breit-Wigner curve of width T, which is extracted from 
data 
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The exact off-shell kinematics and off-shell sub-process 
cross section (at leading order in as) are used. The 
intrinsic-fcr approach is a limiting case of our model at 
r -> 0. 

The model was tested in p| by comparing the triple dif- 
ferential cross section d 3 'a / 'dM 2 dp\dx p of the processes 
pp — > p + fi~ + X and pd — -> /i + pT + X to the data of the 
experiments E866 jl] and E772 0. Here, M denotes the 
invariant mass of the produced lepton pair, xf its Feyn- 
man variable and pt its transverse momentum. Both the 
shape and magnitude of the observed triple differential 
cross sections were well reproduced for reasonable values 
of r w 200 MeV without a need for a if- factor. 

We apply the formalism derived in Q] to predict the 
cross section 
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FIG. 2: Comparison of the predictions of different models: 
(solid) intrinsic-fcT approach; (dash) intrinsic-fcT approach 
with collinear sub-process cross section; (dash-dot) PYTHIA 
generator, xf — 0, D = 0.5 GeV, y/s = 5.5 GeV. The lines 
marked with (1) correspond to M = 1 GeV; (2), M = 2 GeV; 
(3), M = 3 GeV. 



of the process pp — > l + l X in the kincmatical regime 
planned to be probed by PANDA. We use 



da 



dxpdp\dM 2 



flavors 



da 



dxFdp\dM 2 
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where \fs is the centre of mass energy of the pp, the par- 
ton level cross section do is taken from Q], the spectral 
function A is given in (J2J , the unintegrated parton distri- 
bution / is a product of and the parton distribution 
functions [Tl) , the quantities referring to anti-quarks are 
denoted with bars. We assume that T = T and D = D. 

The two parameters of the model (D and T) should 
be extrapolated to the values of M = (1 — 4) GeV and 
Vi = 5.5 GeV relevant for PANDA. In we have per- 
formed fits to the Fermilab experiments E866 and 
E772 9] in different bins of M. From that, we estimate 
the parameter D w 0.6±0.18 GeV. The large uncertainty 
of the extrapolation is reflected in our estimate of the pa- 
rameter error. The model parameter D should be under- 
stood as representing the summed effect of the transverse 
motion of partons inside the nucleon and of the perturba- 
tive corrections 0: D 2 = D 2 ntl +D^ eTt . The transverse 
momentum coming from higher orders of collinear per- 
turbation theory shows at constant r = M 2 / s a linear 
dependence on s in addition to the logarithmic depen- 
dence on M 2 0: (p T ) peTt = s a s (M 2 ) f (r). On the 
other hand, the non-perturbative -Di n t r does not show a 
strong dependence on s. An analysis using the event gen- 
erator PYTHIA [|| has shown that the perturbative cor- 
rections constitute only 0.1% of the pT-spectrum of Drell- 
Yan pairs with pt ^ at PANDA kinematics. Thus, we 
set D « D intr . 

Fig. |21 presents a comparison of the predictions ob- 
tained with three different methods: 1) our model at 
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FIG. 3: Prediction for the pr-distribution of Drell-Yan lepton pairs at PANDA. Intrinsic transverse momentum dispersion 
D = 0.6 GeV. The solid line is the result of calculations in the intrinsic-fc-r approach (width T = 0). The other curves are 
generated with V in the range that we have determined by fitting existing Drell-Yan data. Note that the scale changes from 
nb in the left figure to pb in the right one. -^5 = 5.5 GeV, if = 0.1 in all plots. 
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FIG. 4: Same as Fig. |3] but for higher masses of Drell-Yan pairs. 



r = 0, D = 0.5 GeV, which is equivalent to the intrinsic- 
kx approach (solid lines); 2) the widely used simplified 
version of the intrinsic-fcT approach, in which the depen- 
dence of a on kr is neglected (dash line); 3) PYTHIA, 
taking into account QED and QCD initial state radi- 
ation and intrinsic-fcT (dash-dotted line). The models 
agree in the overall strength of the cross section within 
a factor of 3. The demand of positive mass of the rem- 
nant Mx determines the maximum px- In PYTHIA, 
however, a stronger constraint is implemented: Mx > 
1.6 GeV, where m qq is a diquark mass. We 
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also used 1.6 GeV as the lower bound for the remnant 
mass in the other models for consistency. The dilepton 
mass bin width in the simulation was set to 100 MeV. 
From Fig. one observes that the calculations within 
the simplified intrinsic approach (dashed line) and within 
PYTHIA (dash-dotted line) give similar results. Indeed, 
the intrinsic-fcr approach with an approximate collinear 
sub-process cross section is implemented in PYTHIA. 
The difference between the two curves (dash and dash- 
dot) at high pt values are due to a more detailed mod- 
elling of the remnant in the event generator. Note the 
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FIG. 5: Variation of the cross section with D at fixed F. 
s = 32 GeV 2 , = 0.1, M = 4 GeV. 



qualitative difference of the cross section at the Drell-Yan 
pair mass M = 1 GeV in the intrinsic-fey approach (solid 
lines). In contrast to the higher mass bins, the peak of the 
PT-distribution for M = 1 GeV is not at zero in our model 
(solid line in Fig. [21 • This behavior appears at M < 2D, 
i.e. in the distribution of the low virtuality photons, 
which can be produced by the partonic transverse mo- 
tion alone. It is worthwhile to stress that this drastic 
change in the py-dependence of the cross section takes 
place for all values of T. An experimental verification of 
this effect would constitute a direct test for the transverse 
momentum distribution of partons. Note, however, that 
this effect appears only if the fey-dependence of the par- 
tonic cross section is not neglected (this is the difference 
between the solid and dashed lines in FigEJ) ■ 

Fig. |3 and 0] show several theoretical curves for the 
cross section © with T = (100 - 250) MeV, which is 
the range determined from fitting E866 data in 1]. The 
results of our calculations in the intrinsic-fey approach 
(r = 0) are plotted for comparison (solid lines). The 
evolution of the spectral function with the hard scale is 
unknown. It cannot be directly related to the evolution 



of the fey-distribution, because the quark off-shcllness de- 
pends also on k~. One can see that the effect of the P 
variation on the cross section is of the order of 10%. On 
the other hand, the variation of the parameter D within 
theoretical uncertainty at fixed T (as presented in Fig.lSJ) 
also leads to considerable changes of the cross section. 
As was shown in .1], though, one can do a double fit and 
extract both parameters from the same data set. 

The measurement of this cross section at y/s as low 
as 5.5 GeV will provide essential information on QCD in 
the regime where effects beyond leading order and leading 
twist are expected to be large [14|. The high amplitude 
of the predicted cross sections (up to 10 nb) indicates 
that PANDA at the design luminosity has a potential to 
measure the triple differential unpolarized cross section of 
pp — ► +X with high statistics and an unprecedented 
accuracy. 

The transverse momentum spectrum of Drell-Yan pairs 
at low yfs is generated predominantly by the non- 
perturbative primordial intrinsic transverse momentum 
of the partons. On the other hand, the distribution of 
the intrinsic fey in this region is poorly known. PANDA 
data will be a valuable input that should allow one to 
pin down the quark transverse momentum distribution 
in the proton. Again, the M = 1 GeV mass bin shows 
the mentioned shift of the peak towards higher py. 

The presented results also suggest that one can use 
the future PANDA data to gain information on the spec- 
tral function of partons bound in the nucleon. Indeed, 
as we see from Fig. [21 and 01 an experimental accuracy 
of 20-30% would be enough to answer the question of 
whether the cross section can be described by a model 
with on-shell quarks. The results of the earlier calcula- 
tions at higher y/s together with the PANDA data can be 
used also to extract the dependence of the quark spectral 
function in the proton on the hard scales M and y/s. For 
this purpose, an experimental accuracy of at least 10% 
is needed, so that one can reliably extract the parameter 
r in different mass bins. Should the accuracy be even 
better, the data can be used to investigate the details of 
the quark spectral function shape. 
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